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Thermoacoustic energy conversion is based on the Stirling cycle. In their most basic forms, thermoacoustic 
devices are comprised of two heat exchangers, a porous medium, both placed inside a resonator. Work is 
created through the interaction of strong sound waves with the porous medium that is subject to external 
heating. This work explores the effect of resonator curvature on the thermoacoustic effect. A CFD analysis of a 
whole thermoacoustic engine was developed and the influence of a curved resonator on the thermoacoustic 
effect is discussed. The variation of pressure amplitude and operating frequency serves as metrics in this 
investigation. It was found that the introduction of curvature affects the pressure amplitude achieved. 
Severely curved resonators also exhibited a variation in operating frequency. 

© 2009 Elsevier Ltd. All rights reserved. 


1. Introduction 

The thermodynamic cycle occurring in thermoacoustic heat 
engines (TAEs) as well as thermoacoustic refrigerators (TARs) is the 
Stirling cycle, which was developed in 1816 by Robert Stirling [1 ]. The 
original mechanical Stirling engine utilized two pistons and a 
regenerative heat exchanger [2]. Fig. 1 shows such an engine. 

Over the course of one cycle, the working gas is compressed, and it 
then gives off heat to the heat sink, thus maintaining a constant 
temperature. Afterwards, the gas is heated at constant volume by the 
regenerator and then is heated further at the heat source. This heat 
supply occurs while the gas is allowed to expand and driving the 
power piston, again at constant temperature. After expansion, the gas 
is displaced to the heat sink, while cooling off at constant volume by 
depositing heat to the regenerator, which stores heat between cycle 
segments [2]. The transition to thermoacoustic technology occurred 
when Ceperley recognized that sound waves could replace the pistons 
for gas compression and displacement [3]. 

1.1. Thermoacoustic engines 

The simplest thermoacoustic engines (TAEs) operate with a standing 
acoustic wave in a resonance tube. This wave causes pressure variations 
as well as gas displacement. In a quarter-wavelength resonator, com¬ 
prised of a tube with one open end and one closed end, the pressure 
anti-node of the wave is located at the closed end. The velocity anti-node 
is located at the open end. When a gas is subject to such a standing wave 
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and it interacts with a solid wall that is subject to a (externally imposed) 
temperature gradient, pressure disturbances can be amplified to reach 
strong amplitudes. In practice, this interaction occurs within a porous 
regenerative unit commonly referred to as the “stack” in standing wave 
devices. Ceramic monolith structures are often used as stacks because 
they offer a large surface-to-volume ratio and exhibit desirable hy¬ 
draulic performance. In order for amplification to occur, this temper¬ 
ature gradient within the porous medium must be larger than the 
critical temperature gradient. This critical temperature gradient is 
determined by the temperature changes the gas would undergo if it 
were under the influence of a sound wave in adiabatic conditions. The 
expression for this critical temperature gradient VT crit = c ° Pl was 

PmCpUi 

derived by Swift [4]. It depends on the operating frequency to, the first- 
order pressure and velocity in the standing wave pi and ip, as well as 
the mean gas density p m and specific heat c p . 

In addition to the simple standing wave engine described above is 
the traveling wave engine. The difference between the standing wave 
and traveling wave engine lies in the phasing between velocity and 
pressure. In standing wave engines, the phasing is such that we need to 
artificially delay the transfer of heat by utilizing a porous structure with 
large channels [5] such as the ceramic monolith mentioned previously. 
In the traveling wave engine, the phasing is such that the heating can 
inherently occur after the compression and the cooling can occur after 
the expansion. Here, the flow channels can be designed much smaller 
[5]. Consequently, the temperature of the gas is almost always the 
same as the wall temperature, resulting in heat transfer over very small 
temperature differences. This process produces inherently less entropy, 
and is thus more efficient than the standing wave engine [6,7]. In 
summary, thermoacoustic engines create strong pressure oscillations 
inside a resonator. These oscillations can be used during a reversed 
Stirling cycle to provide cooling. 
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Nomenclature 

c Specific heat 

/ Frequency 

p Pressure 

T Temperature 

u Velocity 


Greek symbols 
A Change 

A Wavelength 

co Angular frequency 

p Density 


Subscripts 

1 First-order 

crit Critical 

m Mean, time averaged 

p Constant pressure 


1.2. Thermoacoustic refrigeration 

If strong sound waves interact with a porous stack their amplitude is 
attenuated, which withdraws thermal energy from the stack. Thermal 
energy is pumped from the stack’s cold end (of temperature below 
ambient) to the exhaust end. Refrigerators and chillers driven by TAEs 
are a reality today; however they are limited to few, specialized uses, for 
example in gas liquefaction. The National Institute of Standards and 
Technology (NIST) in collaboration with Radebaugh built a TAR with 
5 W of cooling power at 120 I< and a low temperature at no load of 901< 
[8]. The main benefit is the lack of moving parts (such as seals), thus 
reducing maintenance costs. It is noteworthy that TARs can achieve 
these low temperatures in a single stage, whereas VCRs can only achieve 
approximately 230 K in a single stage [8]. Cryogenic cooling is not the 
only application for TARs. A practical example of this design was given 
by Poese with a freezer for ice cream storage. This small scale chiller 
featured an annular space around the regenerative unit to achieve 
traveling wave phasing [9]. 

The main reason, why TARs have not been implemented on a large 
scale is based on their inherently low coefficient of performance 
(COP), which is the ratio of input energy over cooling capacity. Note 
that it is not an efficiency in the classical sense and thus not bounded 



by 1. One benchmark of cooling illustrates a spread of COP of tradi¬ 
tional refrigeration between 2 and 6, varying with different working 
fluids and compressor efficiencies [10]. On the other hand, TARs 
are still operating at COPs of less than one [11]. One reason for this 
limitation is that the compression of the gas in TARs is more energy 
intensive than the compression of the liquid working fluid in VC 
refrigerators. Another limitation of the wide spread implementation is 
the physical size of thermoacoustic chillers. If the resonator could be 
coiled as opposed to being designed as a straight tube, the footprint 
of TAEs and TARs could be reduced and thus be advanced towards 
broader implementation. The effect of resonator curvature on the 
performance of TAEs and TARs is not understood. This work provides 
insights into the effect of curvature in general as well as the severity of 
curvature on the thermoacoustic effect and operating frequency. 

2. Development of the CFD model 

In order to investigate different resonator shapes quickly, we 
developed a computational fluid dynamics (CFD) simulation of a stand¬ 
ing wave engine in Fluent. It generally follows the ideas introduced by 
Nijeholt et al. [12] and Hantschk and Vortmeyer [13]. 

2.1. Previous numerical simulations in thermoacoustics 

Nijeholt used ANSYS CFX to simulate a traveling wave engine 
inside a Flelmholtz resonator, mimicking the design built by Bastyr 
and Keolian [14]. They used a very coarse grid, and a time step that 
was just large enough to resolve the expected oscillations. Flantschk 
and Vortmeyer [13] used Fluent 4.4.4 to simulate a Rijke tube (which 
is similar to a classic thermoacoustic engine, except that the oscil¬ 
lations are created with one heated wire screen rather than a stack). 
In their boundary conditions the gas is given a thermal conductivity 
that is 10 times higher than the realistic value for room temperature. 
This corresponds to a Prandtl number that is 10 times lower than the 
real value for their gas at room temperature (which equals approxi¬ 
mately 0.7). Entezam et al. [15] also utilized a CFD code to simulate 
thermoacoustic oscillations in a Rijke tube. In their case applied to a 
pulse combustor. The most recent previous example of a CFD simu¬ 
lation used to investigate the thermoacoustic effect was provided 
by Zoontjens [16]. This work considered a single “thermoacoustic 
couple”, i.e. the interaction of a single stack plate with an oscillating 
gas, as investigated by Ozoe et al. [17]. Zoontjens’ model, was based 
on previous numerical simulation efforts (without using commercial 
codes) put forth byCao etal. [18],Besnoin and co-worker [19,20] and 
Worlikar and Knio [21 ], Ishikawa and Mee [22], and Marx and Blanc- 
Benon [23-26], and finally Piccolo and Pistone [27]. These efforts 
were used largely to investigate the non-linear behavior of the 
temperature as a result of the interaction of the oscillations with 
the stack wall. Our group advanced the numerical simulation by 
accounting for several thermoacoustic couples through a simulation 
of a whole TAE. 

2.2. CFD simulation of a whole TAE 

Our model is based on a mechanically simple standing wave 
engine with a quarter-wavelength resonator and a porous stack 
located in the vicinity of the closed end. This standing wave engine 
was replicated in Gambit. The grid dimensions are 12 mm wide and 
150 mm long. The stack of the engine is represented by several walls 
spaced 500 pm apart. These walls are given a non-zero thickness in 
order to participate in a heat exchange with the surrounding fluid. The 
grid is built using triangular cells, with a total cell count of approx¬ 
imately 36,000 cells. Triangular cells were chosen because it allowed 
to gradually increase the cell density (i.e. the node spacing on the 
walls by the stack is much smaller than on the closed wall boundary 
on the left) towards the stack region, where high accuracy in the 
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Fig. 2. The CFD mesh used for the TAE model. 


representation of the flow behavior is most important. The mesh is 
shown in Fig. 2. It was found that the relative change in the pressure 
amplitude of the thermoacoustic oscillations from an engine with a 
curved resonator compared with the straight case is independent of 
the used grid. However, the qualitative nature of the results presented 
below must be stressed. 

2.3. Boundary conditions 

The simulation of our thermoacoustic engine was performed in 
Fluent Version 6.3.26. The model includes an open end (a 0 Pa gauge 
pressure “pressure outlet”) and a “pressure inlet” (for the initial dis¬ 
turbance) of a positive gauge pressure of 10 Pa. This creates a pressure 
gradient and a non-zero velocity distribution throughout the domain. 
In the present setup, the simulation uses adiabatic walls at the com¬ 
pliance, resonator and stack. The horizontal walls in the stack section 
are given a constant (with respect to time) temperature with a 
gradient from 700 K down to 300 K over a distance of 10 mm utilizing 
a user defined function (UDF). The heat transfer coefficient bet¬ 
ween the stack walls and the fluid is set to be 50 W/m 2 K, which is 
consistent with calculations performed by Liu and Garrett [28]. The 
equation of state used for the working gas is the ideal gas equation, 
which allows for temperature and pressure dependent variations 
of the gas’ density. Viscous and turbulence effects are accounted 
for using the k-s model. The time advancement is first-order, the 
pressure-velocity coupling uses the PISO scheme, and the discretiza¬ 
tion for all variables (pressure, density, momentum and energy) is 
second-order upwind. For the unsteady solution the pressure inlet 
condition on the compliance side was replaced by another adiabatic 
wall, to mimic the real engine. Although previous examples of CFD 
simulations exist in literature, the presented model is the first one to 
be used as a design tool through the investigation of various resonator 
geometries. 

3. Results 

With these boundary conditions in place, we were able to replicate 
the thermoacoustic effect with Fluent. After transitioning from the 
initial pressure disturbance, the simulation yielded strong sustained 
pressure oscillations as shown below in Fig. 3. The same figure also 
contains results from the first comparison between two different 
instances of resonator curvature. Thus, a model of an entire thermo¬ 
acoustic engine with only heat as an input as well as its application as 
a tool to advance the designs of future thermoacoustic engines has 
been demonstrated. Below, the results from the variation of resonator 
curvature will be further elaborated on. 

3A. Effect of resonator curvature 

Again, Fig. 3 shows the transition of the oscillations from the initial 
disturbance to the “limit cycle” (as referred to by Hantschk and 
Vortmeyer [ 13 ] ) as well as the oscillation’s behavior for both cases, the 
straight resonator and a right-angle bend. To account for curvature, 
we utilized a modified grid that was curved 90° over a distance 
of 80 mm (measured from the open end). Thus, the total length of 
the engine remained constant, in order to ensure equal resonance 
conditions and operating frequencies. The resonator was the only 
section that was modified for this investigation, and the node count on 
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Fig. 3. Effect of a curved resonator on thermoacoustic oscillations, comparison between 
0° and 90°. 


the sides of the resonator was maintained as constant for all cases, 
which ensured a close to identical mesh count between the straight 
and curved geometry. As a result, a different pressure behavior is a 
result of the curvature and not from mesh-specific differences. As a 
result of the introduced curvature, we can see that the total amplitude 
of the pressure oscillations is decreased. A close-up of the fully 
developed region is shown in Fig. 4. It becomes obvious that for the 
case of a gradually curved resonator, the pressure amplitude is 
affected, but the frequency of oscillations remains the same. This is not 
the case if the curvature of the resonator is applied over a shorter 
distance (i.e. the resonator is subject to more severe curvature). 

3.2. Variation of severity of curvature 

Above, it was shown that the introduction of curvature causes a 
decrease of pressure amplitude achieved by the TAE. While the 90° 
bend was initially applied over a distance of 80 mm, we also inves¬ 
tigated the same degree of curvature applied over a length of 40 mm. 
Fig. 5 shows a detailed view of the steady oscillations. We can see that 
a change in the severity changes the frequency of the oscillations. 
While the smooth curvature of 90° applied over 80 mm yields a 
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Fig. 4. Detailed view of the effect of curvature (0° and 90°). 
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Fig. 5. Effect of the severity of resonator curvature on thermoacoustic oscillations (90° 
applied over 80 mm and 40 mm). 


frequency of/=614 Hz as in the previous cases, the same curvature 
applied over 40 mm yields an increased frequency of f— 629 Hz. It 
seems to mimic a shorter resonance tube, thus increasing the fre¬ 
quency. The actual tube length now does not directly correspond to 
the operating frequency. On the other hand, the amplitude of the 
oscillations is not affected as the severity of curvature increases. 

4. Conclusion 

In order to advance the application of thermoacoustic devices, 
their footprint should be decreased. We identified the introduction of 
curvature to the resonator as one approach to achieve this reduction 
of the footprint. To investigate the effect of curvature on the per¬ 
formance of a thermoacoustic engine, we have developed a successful 
simulation of said engine in Fluent. Using a small pressure disturbance 
and a driving temperature gradient across the modeled stack as the 
only energy input, we were able to replicate the thermoacoustic effect 
through the amplification of pressure waves. Our model moved 
beyond the investigation of an individual thermoacoustic couple and 
was utilized to quantify (qualitatively) the influence of resonator 
curvature on the engine’s ability to sustain strong thermoacoustic 
oscillations. As a conclusion, we showed that curvature in the reso¬ 
nator influences both the amplitude and the frequency of the sound 
waves, depending on the severity of curvature. Both effects must be 
considered when thermoacoustic refrigerators are designed, because 
a small change in operating conditions can result in drastic changes in 
performance. The effects are summarized below (pressure amplitude 
p [Pa], frequency/[Hz], respective changes denoted by A [%]): 

• Straight resonator (reference case): p — 8584 and/= 614 

• 90°, 80 mm curved: p = 7698, A p = - 10.3,/= 614, A/= 0 

• 90°, 40 mm curved: p = 7625, Ap = -11.2, /= 629, A/= + 2.4. 
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